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Abstract 

The use of electrostatic force microscopy (EFM) to characterize and manipulate surfaces at the nanoscale usually 
faces the problem of dealing with systems where several parameters are not known. Artificial neural networks 
(ANNs) have demonstrated to be a very useful tool to tackle this type of problems. Here, we show that the use of 
ANNs allows us to quantitatively estimate magnitudes such as the dielectric constant of thin films. To improve thin 
film dielectric constant estimations in EFM, we first increase the accuracy of numerical simulations by replacing the 
standard minimization technique by a method based on ANN learning algorithms. Second, we use the improved 
numerical results to build a complete training set for a new ANN. The results obtained by the ANN suggest that 
accurate values for the thin film dielectric constant can only be estimated if the thin film thickness and sample 
dielectric constant are known. 
PACS: 07.79.Lh; 07.05.Mh; 61.46.Fg. 
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Background 

When electrostatic force microscopy (EFM) [1-6] is 
working at the nanoscale, several interacting parameters 
have a strong influence in the signal [7]. Since the elec- 
trostatic force is a long-range interaction, macroscopic 
parameters such as the shape of the tip or the sample 
thickness can strongly modify the electrostatic inter- 
action [8,9]. However, in many experimental situations, 
it is not possible to obtain accurate values for all of these 
parameters, and it is very difficult to achieve quantitative 
experimental results [10]. Previous results [11] have 
shown that artificial neural networks (ANNs) [12] are a 
useful tool to characterize dielectric samples in highly 
undetermined EFM systems. Using known force vs. dis- 
tance curves as inputs for their training, ANNs have 
been able to estimate the dielectric constant of a semi- 
infinite sample in a system where the tip radius and 
shape were not known. 

In this paper, we demonstrate that ANNs can be used to 
improve the accuracy of numerical simulations in EFM 
and to quantitatively estimate the thin film dielectric 
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constant from vertical force curves. First, we compare 
standard minimization and ANN techniques, demonstrat- 
ing that ANN techniques provide a better control of the 
final result of the simulation. The improved numerical 
results are also used to create a complete training set of 
an ANN that estimates the dielectric constant of a thin 
film placed over a dielectric sample. 

As it has been shown before [11], ANNs are able to es- 
timate physical magnitudes in highly undetermined sys- 
tems. In this article, we train an ANN with a complete 
thin film sample to study the necessity of knowing the 
geometry of the sample in the estimations of the thin 
film dielectric constant. Although the influence of the 
thin film thickness is much larger than that of the sub- 
strate dielectric constant, we demonstrate that accurate 
values of the thin film dielectric constant can only be 
obtained when both magnitudes are known. 



Methods 

Artificial neural network formalism for the calculation of 
electric fields 

To briefly illustrate how ANNs can be related to the 
problem of estimating unknown parameters in EFM set- 
ups, let us focus in the scheme shown in Figure la. Here, 
we have a set of metallic objects that are connected to a 
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Figure 1 Two metallic objects with different applied voltages and electrostatic force microscope tip and sample, (a) Schematic 

representation of a system of two metallic objects with different applied voltages V } and V 2 . Its equivalent ANN is also shown as an inset, (b) 

Scheme of an electrostatic force microscope tip and sample. The tip surface has been divided in three regions with a finite number of points (N 1; 

N 2 , and N 3 ). The tip and sample are characterized by eight parameters: the tip sample distance D, the tip apex radius R, the tip length L, the cone 

half-angle 6, the radius R 2l the thin film thickness h h the thin film dielectric constant e h and the sample dielectric constant e 2 
\ J 



battery that provides a constant electric potential. The 
calculation of electric magnitudes such as the electro- 
static potential or the force between these elements is, 
in general, very difficult, and only a few specific geom- 
etries can be analytically calculated [13]. To solve elec- 
trostatic problems with arbitrary geometries, an 
algorithm called the generalized image charge method 
[14,15] (GICM) has been developed. The GICM replaces 



the surface charge density by a set of charges inside the 
metallic objects. The value, position, and number of 
charges are obtained after a standard least-squares 
minimization (LSM) routine for the electrostatic poten- 
tial at the metallic surfaces. An alternative to the LSM is 
to use the ANN formalism by considering the value of 
the charges q x as the weights w x and the potential at the 
metallic surfaces VJ as the expected output values y } (see 




Figure 2 Potential distribution around an EFM tip. 3D representation of the potential distribution around an EFM tip {L = 10 R, 0 = 20°) over a 
dielectric sample (D = 0.1 R, e= 10) as a function of the ANN iterations in the learning process N ]t . 2D equipotential distribution in the middle of 
the image corresponds to the results obtained by the standard LSM 
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Table 1 Coefficients obtained by the ANN and LSM 
algorithms for an EFM system 


Coefficients 


ANN 


LSM 


Q\ 






ql 


-0.26 


417.40 


q2 


1.75 


-2.31 


q3 


-1.54 


2.43 


q4 


0.98 


-0.82 


L } 






L1 


-0.032 


-12,582.71 


L2 


-0.15 


2,194.47 


L3 


-0.30 


-157.42 


L4 


0.90 


5.51 


L5 


0.41 


0.015 


L6 


0.13 


0.48 


L7 


1.32 


-0.39 


L8 


-7.31 


26.26 


L9 


10.16 


-445.25 


L10 


11.51 


3,126.11 


L11 


7.75 


-9,298.52 


L12 


4.42 


9,311.48 



L = 1 0 R, 6 = 20°, D = 0.1 R, and e = 1 0. q t and L t correspond to punctual charges 
and charged segments inside the tip, respectively. 



inset of Figure la). The input patterns Xy play the role of 
the Green functions G^jtyOi), where is the distance be- 
tween the i-charge element and the j -surface point, o- 
represents the geometrical parameters that may be used 
to adequately calculate the electrostatic potential gener- 
ated by the i-charge element (for example, if the element 



is a charged line, o x would represent the length of the 
line). Following this formalism, the electrostatic potential 
Vj can be expressed as 

N c N c 
i=l i=l 

where N c is the number of charged elements q x inside 
the tip. The most right-hand-side term in Equation 1 
represents the electrostatic potential in the notation of a 
single-output ANN, where x% represents the inputs to 
the output neuron y v and w x are the connection weights 
from the inputs (i = 1,. . .,) to this neuron (see Figure la). 
A neural network learns by example. The task of the 
learning algorithm of the network (i.e., the delta rule, 
backpropagation, etc. [12,16]) is to determine w x (i.e., qi) 
from available x^ data. Previous knowledge can help us 
to decide the best values for Af c (by the selection of the 
number of neurons) and G x (by the selection of input 
patterns). 

To compare both minimization techniques (LSM and 
ANN), we have simulated the EFM shown in Figure lb. 
To illustrate the advantages of using the ANN 
minimization routine, we have calculated the q x coeffi- 
cients for the tip-sample distance described in Figure lb 
with both the ANN and standard LSM routines. We 
have used the winGICM software vl.l [17] which also 
uses ANNs to estimate the number of punctual charges 
and the number of segments (4 and 12, respectively). 
The LSM lowest error (located at # = 0.9876 R, 
z = 0.8896 R) at the tip surface was 0.0019 V 0 , where V 0 
is the voltage applied to the tip. In Figure 2, we show 
the electrostatic potential distribution obtained for 




Figure 3 Scheme of the ANN used to estimate e v Inputs are (a) Fvs. D curves; (b) Fvs. D curves and h y thickness; and (c) F vs. D curves, h } 
thickness, and e 2 
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different numbers of iterations N it in the training 
process. We initialized q x = 0 and fixed the learning rate 
to 0.1. When N it = 100,000, the equipotential distribution 
looks identical than that obtained by the LSM. However, 
the highest error at the tip surface is 0.0076 V 0 (four 
times larger than that from LSM). At this point, it seems 
that LSM is a better minimization technique since it 
gives a lower error and does not use any iterative 
process. However, the q { values obtained by the LSM are 
not adequate for several physical applications. As we can 
see in Table 1, the ANN q x absolute values are much 
smaller than those from LSM. This fact is not important 
when q x do not have any physical meaning. However, in 
our case, q x correspond to the charge inside the tip and 
must be used to calculate electric magnitudes like the 
capacitance or the electrostatic force F (used in the fol- 
lowing section). By using ANN q x values, these magni- 
tudes can be calculated with improved accuracy since 



(a) 

100 n 



the low values of the charges strongly reduce the numer- 
ical noise. In conclusion, the ANN minimization allows 
the user to choose the balance between numerical accur- 
acy and physical meaning of the simulations to adapt 
them to the necessity of the problem. 

Results and discussion 

Thin film dielectric constant estimation 

In this section, we are going to use the GICM force vs. 
tip-sample distance (F vs. D) curves for an EFM tip over 
a thin film to train an ANN that will be able to estimate 
the dielectric constant e 2 of the thin film with thickness 
hi (see Figure lb). The thin film will be placed over a 
semi-infinite dielectric substrate characterized by its di- 
electric constant c 2 . It is worth noting that in realistic 
systems where hi is very small, Ci should be considered 
an effective [6] dielectric constant since several nano- 
scale effects can modify the response of the thin material 
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Figure 4 Test error vs. number of iterations and ANN estimation, (a) Test error vs. number of iterations in the ANN learning process for different 
ANN inputs, (b) ANN estimation of the thin film dielectric constant for the test set for the case where Fvs. D, e 2 , and fi, are included as inputs 
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and change the Ci value. Some examples of physical phe- 
nomena that could affect e x are the roughness of the 
thin film surface, the presence of a water layer over the 
thin film [5], or the finite amount of free charge due to 
the small size of the film. The ANN architecture is 
shown in Figure 3. We used a multilayer perceptron 
with sigmoid activation functions. The input layer is 
composed of 20 neurons for the F vs. D curves that are 
calculated for D = {2.5, 5,. . .,50} nm. Additional neurons 
are added in the cases where e 2 and hi are included as 
input values. We used two hidden layers composed of 
10 neurons with no bias applied. The output layer con- 
tains a single neuron which provides the estimate values 
for Ci. We have considered three different inputs: the F 
vs. D curves, h lt and e 2 . GICM F vs. D curves included 
in the training were calculated for e x = {5, 15, 25,. . .,105}, 
e 2 = {5, 15, 25,... ,105}, and /z 1 = {l, 2,. . .,10}. The ANN 
has been tested with 100 F vs. D curves (not used during 
the training) with randomly selected e 1} e 2 (between 5 
and 105), and hi (between 1 and 10) values. As we can 
see in Figure 4a, although the ANN is able to estimate e x 
when F vs. D and hi (excluding e 2 ) are used as input 
parameters, it gives the best results when all the input 
parameters are included. In Figure 4b, we show the error 
obtained by the ANN when all the inputs are included. 
The error is always smaller than 9% for all the e x values. 

The ANN can be used with realistic experimental 
curves without any previous treatment, which is one of 
the advantages of using this technique [11]. In this case, 
experimental curves with a high error could make the 
ANN give wrong e x estimations. This problem can be 
easily solved by training the ANN with a mixture of ex- 
perimental and numerical F vs. D curves. This strategy 
would make the ANN more robust against experimental 
noise (by the use of experimental curves) and still effect- 
ive on the £i estimations (by the use of a whole set of 
numerical curves). 

Recently, a simple analytical expression has been 
developed that demonstrates that a sample composed by 
a thin film over a dielectric substrate gives the same re- 
sponse as that of a semi-infinite uniform dielectric sam- 
ple [18]. The fact that different combinations of e 1} e 2 , 
and hi can correspond to the same effective dielectric 
constant is in agreement with the results found in 
Figure 4a since including e 2 and hi as input values 
improves the ANN performance in the 6i estimations. 

Conclusions 

We have demonstrated that ANNs can strongly improve 
the efficiency of the characterization of samples by elec- 
trostatic force microscopy. First, we have demonstrated 
that the generalized image charge method can be modi- 
fied to use a neural network minimization algorithm. 
Using this technique, we have increased the accuracy of 



the electrostatic force and capacitance calculations. By 
using electrostatic force simulations, we have been able 
to train an ANN to estimate the dielectric constant of 
thin films. The analysis of the results of the ANN sug- 
gests that the thin film dielectric constant can only be 
obtained when the thin film thickness and the dielectric 
nature of the sample are known. Note that the methods 
explained in this paper can be easily applied to experi- 
mental data by providing this kind of input to the ANN. 
If enough data are available, experimental curves can be 
used for the ANN training alone or together with theor- 
etical curves. 
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D: Force vs. tip-sample distance; GICM: Generalized image charge method; 
LSM: Least-squares minimization. 
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